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We study theoretically the effects of finite volume for nn scattering in order to extract physical 
observables for infinite volume from lattice QCD. We compare three different approaches for nn 
scattering (lowest order Bethe-Salpeter approach, N/D and inverse amplitude methods) with the 
aim to study the effects of the finite size of the box in the potential of the different theories, specially 
the left-hand cut contribution through loops in the crossed t, u— channels. We quantify the error 
made by neglecting these effects in usual extractions of physical observables from lattice QCD 
spectra. We conclude that for nn phase-shifts in the scalar-isoscalar channel up to 800 MeV this 
effect is negligible for box sizes bigger than 2.5m^^ and of the order of 5% at around 1.5 — 2m^^. 
For isospin 2 the finite size effects can reach up to 10% for that energy. We also quantify the error 
made when using the standard Liischer method to extract physical observables from lattice QCD, 
which is widely used in the literature but is an approximation of the one used in the present work. 
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I. INTRODUCTION 

One of the aims in present lattice QCD calculations 
is the determination of the hadron spectrum and many 
efforts are devoted to this task |l|-[l9{. A recent review 
on the different methods used and results can be seen 
in [lO]. Since one evaluates the spectrum for particles 
in a finite box, one must use a link from this spectrum 
to the physical one in infinite space. Sometimes, when 
it rarely happens, an energy of the box rather indepen- 
dent of the volume is taken as a proof that this is the 
energy of a state in the infinite volume space. In other 
works the "avoided level crossing", with lines of spectra 
that get close to each other and then separate, is usu- 
ally taken as a signal of a resonance, but this criteria has 
been shown insufficient for resonances with a large width 
[2l| - [23| . A more accurate method consists on the use of 
Liischer's approach, but this works for resonances with 
only one decay channel. The method allows to repro- 
duce the phase-shifts for the particles of this decay chan- 
nel starting from the discrete energy levels in the box 
[13, UHl ■ This method has been recently simplified and 
improved in by keeping the full relativistic two-body 
propagator (Liischer's approach keeps the imaginary part 
of this propagator exactly but makes approximations on 
the real part). The work of [1^ also extends the method 
to two or more coupled channels. The extension to cou- 
pled channels has also been worked out in [26l - [2^ . The 
work of [23} presents an independent method, which is 
rather practical, and has been tested and proved to work 
in realistic cases of likely lattice results. The method has 
been extended in [2^ to obtain finite volume results from 
the Jiilich model for the meson-baryon interaction and in 
[30| to study the interaction of the DK and ijDs system 
where the Z?s*o(2317) resonance is dynamically gener- 
ated from the interaction of these particles. The case of 
the K resonance in the Kw channel is also addressed in 
[3l| following the approach of Ref . [2^ . It has also been 



extended to the case of interaction of unstable particles 
in [3^, to the study of the DN interaction [33)], the tttt 
interaction in the p channel [s^l and to find strategies 
to determine the two A(1405) states from lattice results 

In Ref. jl^ the problem of getting phase-shifts and 
resonances from lattice QCD results ("inverse problem") 
using two coupled channels was addressed. Special at- 
tention was given to the evaluation of errors and the 
precision needed on the lattice QCD calculations to ob- 
tain phase-shifts and resonance properties with a desired 
accuracy. The derivation of the basic formula of |23| is 
done using the method of the chiral unitary approach [36| 
to obtain the scattering matrix from a potential. This 
method uses a dispersion relation for the inverse of the 
amplitude, taking the imaginary part of in th e physi - 
cal region and using unitarity in coupled channels |37l . l38l |. 
The method does not integrate explicitly over the left- 
hand cut singularity. Nevertheless, the latter might lead 
to interesting problems in finite volume calculations be- 
cause in field theory, loops in the t— or w— channel that 
contribute to crossed cuts, are volume dependent. There 
is no problem to incorporate these extra terms into the 
chiral unitary approach by putting them properly in the 
interaction kernel of the Bethe Salpeter equation or N/D 
method l37l . 13911 . or using the inverse amplitude method 
(lAM) j4Cll. I4H. However, the method of [111 to ana- 
lyze lattice spectra and obtain phase-shifts explicitly re- 
lies upon having a kernel in the Bethe Salpeter equation 
which is volume independent. The same handicap oc- 
curs in the use of the standard Liischer approach, where 
contributions from possible volume dependence in the 
potential are shown to be "exponentially suppressed" in 
the box volume. Yet, there is no way, unless one knows 
precisely the source of the volume dependent terms, to 
estimate these effects and determine for which volumes 
the "exponentially suppressed" corrections have become 
smaller than a desired quantity. This is however an im- 
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portant information in realistic calculations. The pur- 
pose of the present paper is to address this problem in a 
practical case, the scattering of pious in s-wave. For that 
we determine the strength of these volume dependent 
terms as a function of the size of the box and the impact 
of these effects in the determination of the phase-shifts 
in the infinite volume case. 

The contents of the paper are as follows. After this 
introduction, we summarize in Sec. |TT] the three models 
used to evaluate tttt scattering in the infinite and finite 
volume case. We then follow by studying the dependence 
on the lattice size of the box L of the resulting phase 
shifts in Sec. IIIII Conclusions are collected in Sec. IIVI 



II. THE TTTT SCATTERING IN THE FINITE BOX 

In this section we explain the three models that we 
are going to consider in the present work to evaluate the 
TTTT scattering within the chiral unitary approach: lowest 
order Bethe-Salpeter (BS), N/D and Inverse amplitude 
method (lAM). The latter two provide contributions to 
the left-hand cut of the scattering amplitude while the 
BS does not. After summarizing the models for the infi- 
nite volume, we explain for each of them how to evaluate 
the scattering in a box of finite size L. We study the 
scalar channel up to total energies of about 800 MeV for 
both isospin (/) and 2. The isoscalar case is relevant 
for the lattice QCD studies of a (or /o(600) meson 
resonance, while for the isotensor case the left-hand cut 
is more relevant (see below). Up to those energies the 
KK and T]r] channels in the 1 = case are negligible, 
hence, we deal here only with the tttt channel. 



A. Lowest order Bethe-Salpeter approach 

In the chiral unitary approach the scattering matrix 
can be given by the Bethe-Salpeter equation in its fac- 
torized form 14311 



T= [1-VG]-W= [V-^~G] 



(1) 



where V is the tttt potential, = — -p(s — ^) for 7 = 
and V = 2j^{s — 2m?) for / = 2, which are obtained 
from the lowest order chiral Lagrangians [4^, with m 
the pion mass and /^r = 92.4 MeV. In Eq. H]) G is the 
loop function of two meson propagators, which is defined 
as 



G 



(2^ (P^ 



1 



1 



•m? + ie 



(2) 



with P the four-momentum of the global meson-meson 
system. Note that Eq. ([T]) only has right-hand cut, unlike 
the other two approaches discussed in the next subsec- 
tions. 



The loop function in Eq. ^ can be regularized ei- 
ther with dimensional regularization or with a three- 
momentum cutoff. The connection between both meth- 
ods was shown in Refs. [ssl, |4l[. In dimensional regu- 
larizatiorQ the integral of Eg. pt, G^, is evaluated and 
gives for the tttt system [s^, llBj 



G^iE) 



(47r)2 



a[fi) +log — +a log ——J 



, (3) 



where a = y 1 — s = E'^, with E the energy of the 
system in the center of mass frame, /x is a renormaliza- 
tion scale and a{fi) is a subtraction constant (note that 
only the combination a{p) — log/i^ is the relevant degree 
of freedom, that is, there is only one independent param- 
eter). 

The loop function G can also be regularized with a 
three momentum cutoff Pmax and, after the integration 
is performed [i^, it results 



G(s) = 



d^p 1 



|pl<Pn 



(27r)3w(p)s-4w(p)2+ie 



(4) 



Let us now address the modifications in order to eval- 
uate the TTTT scattering in a finite box following the pro- 
cedure explained in Ref. [23j . The main difference with 
respect to the infinite volume case is that instead of in- 
tegrating over the energy states of the continuum with p 
being a continuous variable as in Eq. (jj]), one must sum 
over the discrete momenta allowed in a finite box of side 
h with periodic boundary conditions. We then have to 
replace G by G, where 



, |pl<Pmi 



1 



1 



io{p) s — 4:Ll){p}^ 



p = —n , n eZ'^ 

Lj 



(5) 



For the sake of comparison with the other models con- 
sidered in the present work, where dimensional regular- 
ization is always done, we use the procedure of [30| in 
order to write the finite volume loop function G in terms 



^ In our context we refer to the G function given in Eq. (jsjl as 
calculated in "dimensional regularization". Of course, with the 
latter procedure the results is infinite. The infinite is removed by 
the subtraction constant a(/i). A more accurate formulation can 
be given in terms of dispersion relations, the interested reader 
on this point can consult Refs. [stI. |3^ . though the final result is 
the same. 
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of the infinite volume one evaluated in dimensional 
regularization: 



G = G° + lim 



L3 



D<p„ 



(27r)3 



where I{p, s) is the integrand of Eq. Q 



U!(j>) s — 4a;(j?)2 



(6) 



(7) 



Note that G of Eq. depends on the subtraction 
constant a instead of the three- momentum cutoff p,nax- 
In the box the scattering matrix reads 



T 



1 



V- 



G 



(8) 



The eigenenergies of the box correspond to energies 
that produce poles in the T matrix, which corresponds 
to the condition G{E) — V~^{E). Therefore for those 
values of the energies, the T matrix for infinite volume 
can be obtained by 



TiE) = {V-\E) - G{E))-' = [GiE) GiE))' 

The amplitude is related to the phase-shifts by 

SnE 1 



T{E) 



(9) 
(10) 



p cot S — i ' 

where P = -f ^1 — is the CM momentum. 

Eq. © is nothing but Liischer formula 0, HI] except 
that, as shown in Ref. [1^, Eq. © keeps all the terms of 
the relativistic two-body propagator, while Liischer's ap- 
proach neglects terms in Re /(p) which are exponentially 
suppressed in the physical region, but can become sizable 
below threshold, or in other cases when small volumes are 
used or large energies are involved. 



B. The I AM approach 

The next approach considered is the elastic Inverse 
Amplitude Method (lAM) which we briefly review 
in this section and describe how to extend it to consider 
scattering in a finite box. 

The elastic lAM makes use of elastic unitarity and 
Chiral Perturbation Theory (ChPT) ^ to evaluate a 
dispersion relation for the inverse of the tttt scattering 
partial wave of definite isospin / and angular momentum 
J, T" (in the following we drop the superscript IJ to 
simplify notation). The advantage of using the inverse 



of a partial wave stems from the fact that its imaginary 
part is fixed by unitarity. 



ImT = iTp 



ImT" 



Thus, the right-hand cut integral can be evaluated ex- 
actly in the elastic regime and the obtained partial 
wave satisfies unitarity exactly. The partial wave am- 
plitudes calculated in ChPT cannot satisfy unitarity ex- 
actly since they are obtained in a perturbative expan- 
sion T = T2 + + 0{p^), where T2 = 0{p^) and 
T4 = 0{p'^) are the Leading Order and Next-to-Leading 
Order contributions in the chiral expansion of T, respec- 
tively. However, unitarity is satisfied in a perturbative 
way, 



Imr2 ~ 0, Imr4 



-Tr 



2 ' 



(12) 



These equations allow us to evaluate the dispersion rela- 
tion and obtain a compact form for the partial wave as 
we show below. 

We write then a dispersion relation for an auxiliary 
function F = T2 /T, whose analytic structure consists 
on a right-hand cut (RC) from 4m J to 00, a left-hand 
cut {LC) from —00 to 0, and possible poles coming from 
zeros of T, 

F{s) ^ F{0) + F'{0)s + ^F"{0)s^ 

r , ImF(s') (13) 
Jrc s'^(s - s) 

where we have performed three subtractions to ensure 
convergence. In the above equation LC{F) stands for 
the integral over the left-hand cut, and PC stands for 
possible poles contributions, which are present in the 
scalar waves due to the Adler zeros. Using Eqs. pT|) and 
(fT^ we can evaluate exactly in the RG integral ImF — 
— ImT4, and obtain for the right-hand cut RC{F) — 
—RG{T4). The subtraction constants can be evaluated 
with ChPT since they only involve amplitudes or their 
derivatives evaluated at s = 0, F{Qi) ~ r2(0) - ^4(0), 
F'(0) ~ T^(0) - T^(0), F"(Q) ~ -r^'(O). The left-hand 
cut can be considered to be dominated by its low en- 
ergy part, since we have three subtractions, and it is also 
dumped by an extra l/(s' — s) when considering physical 
values of s. Then, we evaluate it using ChPT to obtain 
LG{F) ~ ~LG{T4). The pole contribution is formally 
0{p^) and we neglect it (this causes some technical prob- 
lems in the subthreshold region around the Adler zeros 
which can be easily solved, but they do not affect the de- 
scription of scattering or resonances, for details see [46j). 
Taking into account all the above considerations we ar- 
rive at 



RC{T4) - LC{Ti) - T2[s) - Ti{s) 



(14) 



where in the last step we have taken into account that 
T2{s) is just a first order polynomial in s so that 12(5) = 
T2(0)+T2(0)s, and that the remaining piece in the middle 
member of Eq. is a dispersion relation for —T4{s). 
Then one obtains the simple lAM formula, 



T 



lAM 



rp2 



T2-T4 



(15) 



This formula can be systematically extended to higher 
orders by evaluating the subtraction constants and the 
left-hand cut in the dispersion relation to higher orders. 
Note that the full one-loop ChPT calculation is used, 
so the lAM partial waves depend on the chiral Low En- 
ergy Constants (LECs) , that absorb the loop divergences 
through their renormalization and depend on a renormal- 
ization scale fi. Of course, this fi dependence is canceled 
out in physical observables. In the case of tttt scattering 
there appear four LECs, denoted Z[(^), i = 1...4. These 
LECs are not fixed from symmetry considerations and 
their value has to be determined from experiment. For 
the lAM calculations here we take the values used in [l^l : 
103^5; = -3.7 ± 0.2, 10% = 5.0 ± 0.4, 10^1^ = 0.8 ± 3.8, 
lO^ll = 6.2 ± 5.7, at /i = 770 MeV, which give a good 
description of phase-shift data. Note that in the present 
work we are not interested in a detailed description of 
scattering data, but on the effects of ignoring the expo- 
nentially suppressed dependence on the box size when 
using Liischer's or the chiral unitary approach to obtain 
the scattering phase-shifts from the energy levels in finite 
volume. 

To evaluate the lAM partial waves in a finite box of size 
L we proceed in a similar way as in subsection III Al The 
only piece we need to change is T4, that receives contri- 
butions from loop diagrams, whose momentum integrals 
should be replaced by discrete sums over the allowed mo- 
menta in the finite box. 

Let us first note that the 0{p'^) contribution to the tttt 
scattering amplitude A4{s,t,u), where the partial wave 
amplitude at the same order (r/'^(s)) is obtained by pro- 
jecting ^4(5, t, u) on isospin / and angular momentum J, 
has the form f44l 



^4(5, t, u) ~ B{s, t, u) + C(s, t, u) , 



(16) 



where C(s, t, u) is a second order polynomial in s, t, and 
u which contains the LECs. On the other hand, B{s, t, u) 
contains the non-analyticities coming from the one-loop 
diagrams with vertices from the 0{p^) Lagrangian, 



+ {t{t - u) - 2mlt + Amlu - 2ml} G(t) 
+ {u{u -i)- 2mlu + Amlt - 2ml] G(u) 



(17) 



with 



G{s) 



167r2 



as + 1 „ 

a{s) - 1 



(18) 





G{s) 



G{t) 



FIG. 1. Pion loops in the s (left) and t (right) channels 
contributing to the 0(j)'^) chiral tttt scattering amplitude 
A4,{s,t,u). The s-channel loop gives rise to the right uni- 
tarity cut, whereas the t-channel (as well as u-channel) loop 
contributes to the left-hand cut in the partial wave. 



where we have defined G(s) = G{s) - G(O)0 G(s) is 
the loop function of two pions, Eq. ^ with s = P^. 
There are three loop function contributions, proportional 
to G(s), G{t) and G{u), each one coming from a pion 
loop in the s, t and w-channels respectively, as schemat- 
ically shown in Fig. [TJ The s-channel loops are respon- 
sible for the right unitarity cut, and contain the most 
important L dependence of the amplitude. This L de- 
pendence coming from the unitarity cut is the one used 
by the Liischer/chiral unitary approach method to ob- 
tain the phase-shift from the energy levels in a finite vol- 
ume. However, the t and u-channel loops, which give rise 
to the left-hand cut when projecting into partial waves, 
give an extra dependence on L (polarization corrections 
in the terminology of Ref. 25*|) that is neglected in the 
Liischer/chiral unitary approach method since it is expo- 
nentially suppressed. 

Then, to obtain the lAM amplitudes in finite volume 
we have to replace 24 (s) in Eq. (|15p with T4(s), which is 
obtained projecting into the corresponding partial wave 
the TTTT scattering amplitude in finite volume A4(s,t,u), 
obtained from Eq. but replacing the loop functions 
in Eq. (jl7p with their finite volume counterparts, Gpv- 
Following again the procedure in the finite volume 
loop functions are obtained from the infinite volume ones 
by 



Gfv{z) = G{z) lim 



Hp, z) 



{2TTf' 

(19) 

with z — s, t or u. Now, the energy levels in the 
box are obtained from the poles in the scattering par- 
tial wave (Uni), or equivalcntly, the zeros of 12(3) — T4(s). 



^ This G relates to the J used in [44| as J = —G, in accordance 
with our normalization, where the amplitudes have opposite sign 
to those in | 44j) . 
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From these energy levels at several values of L one can re- 
obtain the phase-shifts for the infinite volume with the 
Liischer/chiral unitary approach method, and compare 
them with the exact infinite volume result to quantify 
the effect of neglecting the L dependence coming from 
the left-hand cut. 



In order to study the finite volume scattering, the same 
replacements as in the lAM and BS methods must be 
done. In particular, in the kernel V{s) — V{s) no change 
is needed in V2(s), whereas Va[s) is changed to Vii-s), 

Vi{s)^fi{s)-T2{sfG{s). (24) 



C. The N/D method 



The case presented in subsection III Al can be put 
in the mor e g eneral framework of the N/D method 
[33, H H, Sg. The amphtude was denoted by T{s) 
in Eq. ([T]). This master formula is obtained by solving 
algebraically the N/D method [13, [H, [H, HO] , with the 
crossed cuts treated perturbatively, while the right-hand 
cut is resummed exactly. The different chiral orders of 
V{s) = T^2(s)+V4(s)-t-. . . are calculated by matching T{s) 
with the perturbative amplitudes T„(s). In this way, up 
to 



T{s) 



V{s) 



1 - V{s)G{s) 

T2{s)+T^[s) + ... 

V2{s) + Vi{s) + V2{sfG{s) + 



(20) 



where the ellipsis indicates 0{p^) 
the expansion. It results then: 



and higher orders in 



V2{s)^T2{s) , 

Vi{s) = Ti{s)-T2{sfG{s) 



(21) 



The finite piece of the unitarity term in the tttt chiral 
amplitude is given by: 



T^{s) = T2{sfG{s) 



(22) 



with G{s) given in Eq. (|18p. In this way, the kernel 
V{s) = V2{s) + Vi{s) has no unitarity cut because: 



TUs) - T2{sfG{s) = T2{sf{G{.s) - G{s)) 



(23) 



and the cut is cancelled in the r.h.s. of the previous 
equation. The full right-hand cut stems then from the 
denominator 1 — V{s)G{s) in Eq. (jT]). 

In the infinite volume case, the LECs are fixed to the 
experiment, as well as the subtraction constant a. We use 
here the central values of the fit given in [5l|, for which 
the values of the finite and scale independent LECs U 
are h = 0.8 ± 0.9, h = 4.6 ± 0.4, U = 2± 4, U = 
3.9±0.5. In terms of the latter, the so-called renormalized 
LECs, which depend on the renornialization scale, are 
IQHl = -2.8 ± 0.9, IQ% = 2.5 ± 0.8, 10% = 2 ± 6, 
10^1^^ = 3 ± 3, where the renormalization scale is chosen 
as fj, = 770 MeV. The subtraction constant a takes the 
value a — —1.2 ±0.4. We additionally note here that the 
same subtraction constant is used for both channels, as 
required by isospin symmetry [53 |. 



Notice that, in view of Eq. there is no effect in 

the s-channel contributions to the kernel V{s). The vol- 
ume dependence enters then in the kernel through the 
t- and M-channel loop functions, where the replacement 
G{z) -> Gfv{z) in Eq. (fTO]! for z = t,u must be done. 
The s-channel volume-dependence enters then at the de- 
nominator of the amplitude T(s) = V{s)/{1 — V{s)G{s)) 
through the function G'(s), Eq. which gives the 

most important contribution to the aforementioned de- 
pendence, as in the case of the lAM method. 



III. RESULTS 

As already explained, the main aim of the present work 
is to quantify the effect of the dependence of the different 
potentials considered on the size of the box, L. Hence, 
we are going to compare the L dependence of the N/D 
and the lAM method with that of the BS, which kernel 
does not depend on L. 




E (MeV) 

FIG. 2. Isospin 1=0, s-wave, tttt — >■ tttt phase-shifts for the 
three different models considered: solid, dashed and dot- 
dashed lines correspond to lAM, N/D and BS, respectively. 
The experimental data are from Refs. [53l - [57| . 
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First we show in Fig. ^ the results for the tttt phase- 
shifts in s-wave and / = for the three different models 
in infinite volume. The lAM and N/D results (solid and 
dashed lines, respectively) are the fits explained in the 
previous section and the BS (dot-dashed line) is fitted in 
this work to the experimental data [ssl - fsTj shown in the 
figure up to 800 MeV. The lAM and N/D approaches are 
essentially equivalent at low energies but differ slightly 
as the energy increases. Thus the difference between the 
lAM and N/D phase shifts is mainly due to the different 
set of data used in the fit and it also gives an idea of 
the theoretical uncertainty. The BS approach produces a 
curve in between the other two, closer to the N/D at low 
energies and to the I AM at higher energies. In any case, 
the different models are compatible within the experi- 
mental uncertainties. Let us note that what matters for 
the discussions in the present work is not the actual val- 
ues of the phase-shifts at infinite volume but the relative 
change when going to the finite box. 

In Fig. Owe show the energy levels for different values 
of the cubic box size, L, for the different models which 
have been obtained from the zeroes of the scattering am- 
plitudes in the finite box as explained in the previous 
section. The dotted lines represent the free tttt energies 
in the box, while the others lines correspond to lAM, 
N/D and BS as in FigM 



900 




L [m ] 



FIG. 3. The first energy levels as a function of the cubic box 
size L for the three different models considered for 1 = 0. The 
dotted lines indicate the free tttt energies in the box. The rest 
of the lines correspond to lAM, N/D and BS as in Fig. [51 



The differences are very small for the largest values of 



L shown in the plot but are more important for smaller 
values of L, specially between the N/D and lAM meth- 
ods. The BS approach produces a curve in between the 
other two, closer to the N/D. The lAM and BS are more 
similar for larger values of energies as can also be seen in 
the phase shifts. Fig. [21 As an example of small L, we 
note that for L = 1.7m~^ the difference between N/D 
and lAM is about 30 MeV. 

An actual lattice calculation would provide some 
points over analogous trajectories in the E vs. L plots. 
The "inverse problem" is the problem of getting the ac- 
tual scattering amplitudes (and hence by-product mag- 
nitudes like phase-shifts) in the infinite space from data 
produced by lattice QCD consisting of points in plots of 
E vs. L over the energy levels in the box. For points in 
these levels the amplitude in the infinite volume can be 
extracted from the generalization of the Liischer formula, 
as explained in the previous sections. 



T{E) 



1 



G{E) - G{E) 



(25) 



100 ■ 



90 - 



lAM infinite volume 
lAM (level 2) 
N/D infinite volume 
N/D (level 2) 



<u 
-a 



to 




,50 300 350 400 450 500 550 600 650 700 750 800 



E (MeV) 

FIG. 4. Solution of the inverse problem for J = for the 
I AM and N/D methods. The BS result is the same as in the 
infinite volume case and thus we do not show it in the figure. 
We show the results obtained only from level 2 of Fig. [31 since 
the results with levels > 2 almost overlap with the infinite 
volume line. For the meaning of each line consult the inset in 
the figure. 



In Fig. m we show the phase-shifts obtained for the 
different methods implementing the "inverse problem" 
analysis (or "reconstructed" results) with Eq. (|25p and 
from the E vs. L plot. For the BS model the results are 
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independent of the level used for a given E, since the po- 
tential does not depend on L, and they are equal to the 
infinite volume result. Therefore we do not show the BS 
result since it is the same as in Fig. [5] For the I AM and 
N/D methods the results depend on the level chosen for 
a given E since the potentials depend on L as explained 
in the previous sections. Actually, for levels > 2 of Fig. [3] 
the results are almost equal to the infinite volume results 
and hence we do not show them in the figure since they 
would almost overlap with the infinite volume line. This 
is because, as seen in Fig. [31 for the higher energies these 
levels imply large values of L. Indeed, for energies be- 
low 800 MeV this implies values of L higher than about 
3to^^. For the results obtained with level 2, the phase- 
shifts differ in about 5% of the result in the infinite vol- 
ume at the higher energies considered. For E ~ 800 MeV 
this implies L values slightly smaller than 2m~^, as can 
be seen in Fig. [31 It is worth noting that the effect of the 
dependence on L of the models with left-hand cut go in 
the same direction and are of similar size in spite of the 
different models used. This gives us confidence that the 
actual L dependence of the left-hand cut is properly con- 
sidered and the real effect of any realistic model would 
be of the order obtained in the present work. An analysis 
with Eq. ([^ applied to actual lattice results of E versus 
L levels would neglect the possible L dependence of the 
potential and hence the errors from the L dependence of 
the left-hand cut would be of the order of the differences 
shown in the figure. Note also that the L dependence of 
the results are smaller than the initial difference between 
the N/D and lAM themselves and also lower than the 
experimental uncertainties. Therefore, an actual lattice 
calculation should care about this L dependence only if 
it aims at getting errors smaller than the effect obtained 
in the present work. 

In Figs. [3 [21 and [3 we show for the 1 — 2 case the 
same results as in Figs. [H to [H for / = 0. In Fig. [51 we 
see that the I AM and N/D methods provide very similar 
results and compatible with the experimental data while 
the BS approach gets worse phase-shifts. This is because 
in the lAM and N/D the left-hand cut is included per- 
turbatively order by order, unlike the BS model, and in 
this channel the left-hand cut is more relevant. In Fig.[^ 
we show the energy levels in the box for this channel. 

Now both lAM and N/D provide similar results. In 
Fig. [3 we show the solution of the inverse problem for 
the phase-shifts. We see that the N/D method provides 
a higher L dependence for large values of the energies, 
unlike lAM. At 800 MeV the difference is about 10% for 
the N/D and 2% for the lAM. The difference in the phase- 
shifts between the two approaches is large in spite of the 
energy levels being very similar. This is because the en- 
ergy levels are very close to the free case, unlike the 7 = 
case, and then the G function is very steep. This makes 
that small variations in E provide large variations in G. 

In usual inverse problem analysis from actual lattice 
results, it is common to use the Liischer formula [H, [25j 
which, as explained in section FlI Al is an approximation 
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FIG. 5. Isospin 1—2, s-wave, tttv — >■ nn phase-shifts for the 
three different models considered. The experimental data are 
from Refs. [s^. [59l]. See the inset in the figure for the corre- 
spondence between the different lines and the approach used. 



to that used in the present work, Eq. Therefore 
it is worth studying what is the error made in the re- 
constructed phase-shifts if one uses the Liischer equation 
instead of the exact one. In Ref. [1^ it was shown that 
the Liischer method can be reproduced if in Eq. Q one 
substitutes 



by 



Hp, s) 



Hp, s) 



1 



1 



PoN - P' 



(26) 



(27) 



where poN-f^l-^. 

In Fig. [51 we show the effect in the isospin phase-shifts 
of using the pure Liischer method, Eq. (|771) . instead of the 
exact one, Eq. (|^ . (For the isospin 2 case the effect is 
small and thus we do not show any plot.) The difference 
is significant only for phase-shifts extracted from level 
2 of Fig. [31 since the difference is only relevant for small 
values of i. Therefore we only plot results extracted from 
level 2. The difference between the exact method and the 
Liischer one is similar for all the three different models 
for the potential. The size of the difference is similar to 
the one from the L dependence of the potential discussed 
above but goes in the opposite direction. Therefore they 
tend to compensate each other by chance. 
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FIG. 6. The first energy levels as a function of the cubic box 
size L for the three different models considered for 1 = 2. The 
meaning of the different lines is as in Fig. [3l 
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FIG. 7. Solution of the inverse scattering problem with 1 = 2 
for the I AM and N/D methods. The BS result is the same as 
in the infinite volume case and thus we do not show it in the 
figure. We show the results obtained only from level 1 and 2 
of Fig. [6] since the results with levels > 2 almost overlap with 
the infinite volume line. For the meaning of the lines consult 
the inset in the figure. 



IV. SUMMARY 



In this paper we have faced for the first time in the lit- 
erature the problem of the presence of the left-hand cut 
of the TTTT amplitude for the evaluation of phase-shifts 
from lattice QCD results using Liischer's approach. The 
t— and u— channel terms can be taken into account in a 
field theoretical approach by means of the lAM, or N/D 
NLO methods, leading to good reproductions of the scat- 
tering data. Results from lattice QCD should contain all 
the dynamics and, as a consequence, should account for 
these effects too. However, the method to go from the 
discrete energy level in a box from lattice simulations to 
the phase shifts for scattering in the infinite volume case 
requires the use of Liischer's approach, or its improved 
version of [1^, both of which rely upon the existence 
of a volume independent potential. Yet, the terms con- 
tributing to the left-hand cut, containing loops in the t— 
and u— channels, are explicitly volume dependent. In this 
work we have investigated the errors induced by making 
use of [i^ or [2^ in the reproduction of phase-shifts from 
the energy spectra of lattice calculations in the finite box. 
We have found that in the case of tttt scattering in s-wave, 
both for / = and 1 = 2, the effect of the L dependence 
in the potential is smaller than the typical errors from the 
experimental phase-shifts or the differences between the 
three models that we have used, the lAM, N/D NLO and 



BS LO. This is good news for lattice calculations since 
one of the warnings not to go to small values of L was 
the possible L dependence of the potential which in some 
cases, like in the present one, we know that exists. We 
found that it is quite safe to ignore this dependence for 
L > 2.5m~^, and even with values of L around 1.5— 2to~^ 
the errors induced are of the order of 5%. 

On the other hand we have quantified the error made 
by using the pure Liischer formula instead of the exact 
one, Eq. The effect in the phase-shifts of this ap- 

proximation tends to compensate, by chance, the effect 
of neglecting the L dependence in the potential discussed 
so far. 

All these findings, together with the use of the ap- 
proach of [11] that also eliminates L depended terms (ex- 
ponentially suppressed) from the Liischer's approach, can 
encourage the performance of lattice calculations with 
smaller size boxes with the consequent economy in the 
computing time. 
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FIG. 8. Difference between the exact inverse method for- 
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approach corresponding to each line is given in the inset in 
the figure. 
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